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Abstract.  We  have  quantihed  the  relationship  between 
Aerosol  Index  (AI)  measurements  and  plume  height  for 
young  biomass  burning  plumes  using  coincident  Ozone 
Monitoring  Instrument  (OMI)  and  Cloud- Aerosol  Lidar  and 
Infrared  Pathhnder  Satellite  Observations  (CALIPSO)  mea¬ 
surements.  This  linear  relationship  allows  the  determination 
of  high-altitude  plumes  wherever  AI  data  are  available,  and 
it  provides  a  data  set  for  validating  global  hre  plume  heights 
in  chemistry  transport  models.  We  hnd  that  all  plumes  de¬ 
tected  from  June  2006  to  February  2009  with  an  AI  value  >9 
are  located  at  altitudes  higher  than  5  km.  Older  high-altitude 
plumes  have  lower  AI  values  than  young  plumes  at  similar 
altitudes.  We  have  examined  available  AI  data  from  the  OMI 
and  TOMS  instruments  (1978-2009)  and  hnd  that  large  AI 
plumes  occur  more  frequently  over  North  America  than  over 
Australia  or  Russia/Northeast  Asia.  According  to  the  derived 
relationship,  during  this  time  interval,  181  plumes,  in  various 
stages  of  their  evolution,  reached  altitudes  above  8  km. 


1  Introduction 

Knowledge  of  the  injection  and  transport  altitudes  of  biomass 
burning  plumes  is  crucial  for  understanding  the  long  range 
transport  of  particles  and  reactive  gases  that  can  substan¬ 
tially  alter  the  radiation  balance  and  chemistry  of  the  atmo¬ 
sphere  (Cofer  et  ah,  1996;  Waibel  et  ah,  1999;  Forster  et 
ah,  2001;  lost  et  ah,  2004;  Luderer  et  ah,  2006;  Rander- 
son  et  al.,  2006;  Trentmann  et  ah,  2006;  IPCC,  2007;  Jeong 
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and  Hsu,  2008).  Remote  and  in-situ  observations  have  docu¬ 
mented  many  instances  of  biomass  burning  pollutants  in  the 
upper  troposphere  and  lower  stratosphere  (UT/LS)  (Fromm 
and  Servranckx,  2003;  lost  et  al.,  2004;  Rosenfeld  et  ah, 
2007). 

To  date,  biomass  burning  plume  heights  have  typically 
been  determined  in  a  case-by-case  manner.  For  instance,  us¬ 
ing  Polar  Ozone  and  Aerosol  Measurement  III  (POAM  III) 
observations  and  back  trajectories  (Fromm  and  Servranckx, 
2003;  Rosenfeld  et  al.,  2007;  Fromm  et  al.,  2008a,  b),  the 
May  2001  Chisholm  hre  in  Alberta,  Canada  was  determined 
to  have  penetrated  the  tropopause.  In  the  northern  summer 
of  2001  over  200  observations  of  aerosol  layers  were  mea¬ 
sured  by  the  POAM  III  mission,  half  of  which  were  in  the 
stratosphere  (Fromm  and  Servranckx,  2003).  Another  case 
study  involved  detection  of  high  CO  values  2  km  above  the 
tropopause  south  of  Florida  during  an  aircraft  campaign  in 
July  2002  (lost  et  al.,  2004).  Particle  measurements  by  the 
Particle  Analysis  by  Laser  Mass  Spectrometry  (PALMS)  in¬ 
strument  helped  corroborate  the  biomass-burning  signature. 
In  addition,  isentropic  back  trajectories  placed  the  air  mass 
under  a  POAM  III  enhanced  aerosol  feature  over  Canada, 
showing  evidence  of  stratospheric  transport  far  from  the  ob¬ 
servation  location. 

Mazzoni  et  al.  (2007)  used  Multi-angle  Imaging  Spectro- 
Radiometer  (MISR)  and  Moderate  Resolution  Imaging  Spec- 
troradiometer  (MODIS)  data  to  locate  hres  and  their  associ¬ 
ated  smoke  plumes  and  retrieved  the  injection  heights  gen¬ 
erated  by  hre  buoyancy  over  a  four-month  period.  Smoke 
plume  injection  height  is  dehned  as  the  altitude  at  which 
the  smoke  particles  are  injected  to  the  atmosphere  before 
transport  (Kahn  et  al.,  2008).  Mazzoni  et  al.  found  a  to¬ 
tal  of  325  candidate  plumes,  of  which  138  were  identihed  as 
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biomass  burning  plumes,  none  of  which  were  classified  as 
stratospheric.  Kahn  et  al.  (2008)  further  assessed  smoke  in¬ 
jection  height  for  a  similar  region  and  period.  They  noted  that 
at  least  about  10%  of  wildfire  smoke  plumes  reached  the  free 
troposphere  but  remained  lower  than  5  km.  This  low-altitude 
bias  may  be  inherent  in  their  selection  criteria,  as  their  plume 
observations  are  limited  to  the  overpass  time  (11:00-14:00 
local  time  in  North  America,  Val  Martin  et  al.,  2010)  of  the 
Terra  satellite,  when  fires  have  not  yet  reached  their  maxi¬ 
mum  intensity. 

Val  Martin  et  al.  (2010)  extended  the  work  of  Mazzoni 
et  al.  (2007)  by  analyzing  a  5-year  record  of  MISR  smoke 
plume  injection  heights  over  North  America  and  correlat¬ 
ing  them  with  MODIS  fire  radiative  power  measurements. 
That  study  highlighted  a  physical  relationship  which  could 
be  explored  to  estimate  injection  heights  for  low-altitude  fire 
plumes  (lower  than  7  km). 

Recently,  Jeong  and  Hsu  (2008)  presented  a  new  algorithm 
that  predicts  the  height  of  smoke  aerosols  associated  with 
biomass  burning  when  Aerosol  Index  (Al),  single  scattering 
albedo  (SSA),  and  optical  depth  are  known.  They  showed 
that  for  a  specific  biomass  burning  plume  over  the  north¬ 
west  United  States  (15  August  2007),  the  predicted  height 
and  the  Cloud- Aerosol  Lidar  and  Infrared  Pathfinder  Satel¬ 
lite  Observations  (CALIPSO)  measured  height  are  in  rela¬ 
tively  good  agreement  (correlation  coefficient  =  0.86).  Their 
results  showed  that  satellite  data  can  be  exploited  in  studies 
of  larger  scope. 

Several  models  have  been  used  to  assess  the  impact  of 
biomass  burning  pollution  on  the  stratosphere  (e.g.  Luderer 
et  al.,  2006;  Trentmann  et  al.,  2006),  the  UT/LS  (e.g.  Dun¬ 
can,  2007)  and  the  troposphere  (e.g.  Kuester  et  al.,  2005;  Li 
et  al.,  2005).  The  injection  height  of  smoke  plumes  from  for¬ 
est  fires  is  a  key  input  for  aerosol  transport  modeling,  as  the 
height  is  critical  for  determining  the  distance  and  direction 
the  smoke  will  travel  (Westphal  and  Toon,  1991;  Ginoux  et 
al.,  2001;  Colarco  et  al.,  2004).  However,  the  injection  height 
is  a  large  source  of  uncertainty  in  these  models.  Recently,  ef¬ 
forts  have  been  made  to  apply  more  realistic  injection  heights 
in  chemistry  transport  models.  These  physically-based  (Fre¬ 
itas  et  al.,  2007;  Guan  et  al.,  2008)  or  MISR-derived  (Chen  et 
al.,  2009)  injection  heights  have  only  been  validated  in  cer¬ 
tain  regions  and  for  limited  periods.  In  order  to  minimize  un¬ 
certainties  and  limitations,  it  is  important  to  use  a  database  of 
seasonally  and  regionally  diverse  plume  heights,  rather  than 
simply  compare  model  results  to  a  handful  of  case  studies. 

We  have  developed  a  simple  empirical  method  using 
the  Al  as  determined  by  the  Ozone  Monitoring  Instrument 
(OMI)  and  Total  Ozone  Mapping  Spectrometer  (TOMS)  in¬ 
struments  to  identify  biomass  burning  plume  heights.  The 
Al  is  a  useful  tool  for  detecting  UV-absorbing  aerosols  since 
it  is  derived  from  the  measurement  of  the  reflectance  at  two 
wavelengths  in  the  UV  coupled  with  radiative  transfer  mod¬ 
eling  of  the  Rayleigh  atmosphere.  The  Al  is  particularly  sen¬ 
sitive  to  high-altitude  aerosols  and  has  a  multi-decadal  histor¬ 


ical  data  set.  We  have  used  plume  altitudes  measured  from 
the  CALIOP  instrument  to  develop  a  method  for  screening 
the  entire  Al  data  record  for  high  altitude  plumes.  Using 
the  relationship  developed  when  altitude  measurements  are 
available,  we  screened  the  entire  Al  record  and  assembled 
a  climatology  of  high-altitude  biomass  burning  plumes  from 
1978  to  2009.  In  addition  to  providing  a  long-term,  global 
data  set  for  model  validation,  our  data  set  also  reveals  the  pat¬ 
terns  and  numbers  of  high  altitude  biomass  burning  plumes 
in  several  regions  worldwide. 

2  Satellite  data  sets 
2.1  Aerosol  Index  data 

The  Al  is  a  measure  of  how  backscattered  ultraviolet  (UV) 
radiation  from  an  atmosphere  containing  aerosols  differs 
from  that  of  a  pure  molecular  atmosphere.'  Al  is  especially 
sensitive  to  the  presence  of  UV  absorbing  aerosols  such  as 
smoke,  mineral  dust,  and  volcanic  ash.  The  index  value  is 
positive  when  absorbing  aerosols  are  present,  whereas  clouds 
yield  a  nearly  zero  value  of  AL  Torres  et  al.  (1998)  and  Jeong 
and  Hsu  (2008)  showed  that  Al  varies  with  aerosol  layer 
height,  optical  depth  and  single  scattering  albedo.  High  and 
optically-thick  absorbing  aerosol  layers  more  effectively  re¬ 
duce  the  scattered  light  from  the  Rayleigh  scattering  below 
the  aerosol  layer.  In  cloud-free  conditions,  the  highest  and 
thickest  absorbing  aerosols  give  the  largest  Al  values.  Ab¬ 
sorbing  aerosols  can  still  be  detected  if  mixed  with  clouds  or 
located  above  clouds. 

The  TOMS  instruments  onboard  the  Nimbus  7  (1978- 
1993)  and  Earth  Probe  (1996-2005)  satellites  and  the 
OMI  onboard  the  EOS -Aura  (2004-present)  satellite  pro¬ 
vide  a  nearly  continuous  and  long  term  record  of  Al  suit¬ 
able  for  scientific  applications.  The  TOMS  and  OMI  Al 
datasets  are  derived  using  the  same  methodology  and  wave¬ 
length  channel  (Ahmad  et  al.,  2006;  Li  et  al.,  2009). 
The  consistency  of  both  datasets  has  been  demonstrated 
in  a  recent  study  (Li  et  al.,  2009),  offering  a  basis  for 

'  Al  is  defined  as: 

AI(r,iu,z)=100(logio(/360meas/631  meas)-loglo(660  calc(^sfc)/631  calc(.Ssfc))), 

where  I  is  the  radiance  (or  intensity)  measured  by  the  satellite  or 
calculated  for  a  purely  molecular  scattering  atmosphere  and  is 
the  Lambert  equivalent  effective  reflectivity  value.  The  variables 
T,  ft),  and  z  represent  optical  depth,  single  scattering  albedo,  and 
height,  respectively.  This  rather  awkward  definition  is  due  to  the 
origin  of  Al  as  a  residual  from  satellite  measurements  (the  term  was 
called  an  “N-valued  residual”  in  the  original  Herman  et  al.  (1997) 
paper).  In  practice  the  equation  reduces  to  the  negative  of  100  times 
the  log  base  10  of  the  ratio  of  the  measured  radiance  at  a  wavelength 
divided  by  the  radiance  that  would  be  observed  if  the  atmosphere 
had  only  molecular  scattering.  The  use  of  Al  has  been  discussed  in 
a  number  of  papers  (Torres  et  al.,  1998;  Jeong  and  Hsu,  2008). 
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Fig.  1.  Locations  of  AI  plumes  encountered  during  the  CALIPSO  era  (1  June  2006  to  13  February  2009).  An  aerosol  plume  is  defined  as 
containing  at  least  1  pixel  with  AI>7  and  at  least  3  adjacent  pixels  with  AI>4.  Orange  and  red  symbols  represent  the  cases  with  and  without 
coincident  AI  and  CALIPSO  measurements,  respectively. 


combined  studies  (Gautam  et  al.,  2009).  Daily  TOMS  ver¬ 
sion  8  Level  3  (1.25°  x  1°  resolution)  and  OMI  Collection  3 
(1°  X  1°  resolution)  were  obtained  from  the  NASA/GSFC 
TOMS  website  (http://toms.gsfc.nasa.gov/aerosols/aerosols. 
html,  updated  June  2009).  TOMS  AI  data  was  unavailable 
between  May  1993  and  September  1996.  Additionally,  the 
EP-TOMS  AI  data  after  July  2000  were  not  included  in  this 
study  due  to  the  calibration  drift  of  the  sensor  (Kiss  et  al., 
2007). 

Over  bright  clouds,  AI  values  are  enhanced  due  to  the  ef¬ 
fect  of  the  radiative  transfer  interactions  between  the  bright 
underlying  cloud  and  the  absorbing  layer  (Hsu  et  al.,  2004; 
Torres  et  al.,  2007)  and  do  not  represent  actual  differences 
in  aerosol  physical  properties.  To  avoid  the  enhanced  effect 
of  underlying  clouds  on  AI,  we  screened  out  cloudy  pixels 
that  contained  a  cloud  fraction  (for  OMI)  or  a  reflectivity  (for 
TOMS)  of  greater  than  0.2  (personal  communication,  O.  Tor¬ 
res,  2008). 

2.2  CALIPSO  data 

CALIPSO’s  lidar  instrument,  the  Cloud-Aerosol  Lidar  with 
Orthogonal  Polarization  (CALIOP),  provides  high  resolution 
vertical  profiles  of  aerosol  and  cloud  attenuated  backscat- 
ter  signals  at  532  nm  and  1064  nm  (Winker  et  al.,  2007). 
The  footprint  of  the  lidar  beam  has  a  100  m  cross-section 
(Thomason  et  al.,  2007)  with  an  overpass  around  1:30p.m. 
local  time.  CALIPSO  gives  limited  coverage  relative  to  OMI 
(swath  =  2600 km)  (Eckhardt  et  al.,  2008). 

Several  products  derived  from  the  backscatter  signals  were 
used  in  this  study.  They  included  the  version  2  Vertical  Eea- 
ture  Mask  (VEM)  (Vaughan  et  al.,  2004),  the  532  nm  total 
attenuated  backscatter,  the  depolarization  ratio  (the  ratio  of 
the  perpendicular  component  of  the  backscatter  to  the  paral¬ 
lel  component),  and  the  color  ratio  (the  ratio  of  backscatter  at 
1064  nm  to  532  nm).  We  used  the  VEM  product  to  determine 


the  top  heights  of  features  and  then  further  distinguished  the 
aerosol  from  cloud  by  checking  the  532  nm  total  attenuated 
backscatter,  depolarization  ratio,  and  color  ratio.  Details  on 
distinguishing  aerosols  from  clouds  are  described  in  Sect.  3. 

2.3  MODIS  fire  data 

Additionally,  MODIS  fire  count  maps  were  used  to  detect 
the  location  and  time  of  fires  for  the  2004-2009  period.  Eire 
count  data  are  available  online  at  http://firefly.geog.umd.edu/ 
flremap/. 

3  Identification  methodology  for  high-altitude  plumes 

We  first  used  the  AI  data  product  to  locate  large  aerosol 
plumes,  defined  as  containing  at  least  1  pixel  with  AI>7  and 
at  least  2  (TOMS)  or  3  (OMI)  adjacent  pixels  with  AI>4. 
We  found  1479  plumes  over  the  AI  record  period  (1978- 
2009,  see  Sect.  2.1  for  specific  dates).  Erom  that  set,  large 
AI  observations  caused  by  solar  eclipses  and  volcanoes  were 
removed.  Major  desert  regions  which  are  likely  influenced 
by  dust  rather  than  smoke,  such  as  North  Africa  and  the 
Middle  East,  were  excluded.  Because  the  dust  storm  season 
overlaps  with  the  wildfire  season  in  Russia,  Mongolia,  and 
Northeast  China,  we  screened  out  potential  dust  cases  over 
the  downwind  regions  of  the  Gobi  and  Taklamakan  Deserts 
during  spring  by  checking  the  MODIS  Are  product  for  the 
time  period  that  a  plume  was  first  observable  or  by  track¬ 
ing  the  plume’s  geographical  origin.  In  addition,  news  re¬ 
ports  and  published  literature  were  used  as  tools  in  linking 
an  aerosol  plume  to  a  dust  signal,  which  was  then  removed 
from  further  analysis.  We  found  107  AI  plumes  likely  due 
to  smoke  aerosol  during  the  CALIPSO  era  (1  June  2006  to 
13  Eebruary  2009).  The  location  of  these  smoke  plumes  is 
displayed  in  Eig.  1 . 
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Fig.  2.  Flowchart  for  the  determination  of  highest  plume  altitude 
along  the  CALIPSO  track. 


Once  the  smoke  plumes  were  located,  we  determined 
the  plume  altitude  following  the  steps  highlighted  in  the 
flowchart  of  Fig.  2.  We  first  used  the  CALIPSO  level  2 
VFM  data  to  find  the  plume  altitude.  If  the  CALIPSO  day 
track  passed  through  the  Al  plume,  we  found  the  highest 
CALIPSO  signal  height  inside  the  Al  plume  area.  An  in¬ 
spection  of  CALIPSO  images  for  all  plume  cases  shows  that 
plume-top  heights  are  relatively  uniform  within  each  plume 
with  a  maximum  variation  within  2  km.  It  should  be  men¬ 
tioned  that  for  an  optically  dense  smoke  layer  discussed  here, 
the  median  height  is  difficult  to  determine,  because  the  low¬ 
est  altitude  from  which  the  lidar  signal  is  returned  may  lie 
well  above  the  true  base  of  a  dense  plume.  Thus,  our  analysis 
addresses  only  plume-top  heights.  In  cases  where  more  than 
one  CALIPSO  track  intersected  the  plume,  the  track  that  was 
closest  to  the  location  of  the  largest  Al  value  was  used.  The 
highest-altitude  signal  observed  by  CALIPSO  at  the  location 
of  a  large  Al  observation  should  be  related  only  to  absorbing 
aerosol,  because  cloud  produces  a  near-zero  Al  value. 

Because  the  Al  grid  may  contain  sub-pixel  sized  clouds, 
which  can  be  explicitly  detected  by  the  high-resolution 
CALIPSO  measurement,  we  also  checked  the  available  532- 
nm  total  attenuated  backscatter,  CALIPSO  depolarization 
ratio,  and  the  color  ratio  images  for  each  case  to  validate 
that  the  high-altitude  signals  did  represent  smoke  aerosols. 
Unlike  clouds,  thin  smoke  plumes  have  a  weak  532  nm  to¬ 
tal  attenuated  backscatter  signal;  dense  plumes  may  have 
a  532  nm  backscatter  signal  similar  to  clouds.  We  distin¬ 
guished  dense  smoke  from  cloud  by  checking  for  an  obvi¬ 
ous  increase  of  color  ratio  with  layer  penetration  depth  cou¬ 
pled  with  low  depolarization  ratio  (dominated  by  a  value  less 
than  ~0. 1)  (personal  communication,  A.  Omar,  2008).  If 
the  top  signal  was  a  cloud,  we  assigned  the  highest  smoke- 
plume  height  within  the  plume  track  as  the  top  height.  The 


Fig.  3.  Scatter-plot  of  OMI  Al  and  CALIPSO  maximum  plume 
height.  Red  and  blue  symbols  represent  “young”  and  “aged” 
plumes,  respectively.  The  inclined  black  line  (Z=  L2AI  —  5.4)  is 
the  lower  bound  to  the  young  plume  data.  The  black  vertical  line 
indicates  an  Al  threshold  (Al  =  9)  for  high-altitude  (>5  km)  plumes. 
The  red  line  (Z=  1.1  Al  —  1.2)  is  the  best  fit  to  the  young  plumes. 
The  shaded  area  shows  ail  standard  deviation  (1 -sigmas^  1.9  km) 
across  the  best  fit  of  the  young  plumes. 


height  uncertainty  resulting  from  this  manual  assignment 
was  approximately  ±0.5  km.  Using  the  above  method,  we 
found  23  cases  of  coincident  Al  and  CALIPSO  measure¬ 
ments  which  are  plotted  in  Fig.  3  and  also  denoted  by  orange 
symbols  in  Fig.  1 . 

By  running  a  back  trajectory  model  starting  at  each  grid 
point  (1.25°  X  1°)  within  the  23  Al  plumes  and  their  associ¬ 
ated  maximum  CALIPSO  plume  height,  we  linked  each  Al 
plume  to  the  corresponding  fire  source  region  observed  by 
MODIS  and  estimated  the  plume  age.  The  trajectory  model 
employed  here  was  the  Goddard  kinematic  trajectory  model 
(Schoeberl  and  Newman,  1995)  using  National  Centers  for 
Environmental  Prediction  (NCEP)  winds  updated  every  6  h 
on  a  2.5°  x  2.5°  grid.  The  maximum  Al  measurement  time 
(local  1:30  p.m.)  was  used  as  the  starting  time  of  trajec¬ 
tory  model.  Plumes  that  can  be  tracked  back  to  a  fire  region 
within  2  days  were  classified  as  “young”,  while  plumes  that 
have  been  removed  from  their  fire  regions  for  longer  than 
two  days  were  defined  as  “aged”.  We  found  19  young  and 
4  aged  plumes,  which  are  shown  in  Pig.  3  as  red  and  blue 
symbols,  respectively.  Table  1  lists  the  date,  maximum  Al 
and  top  height  of  each  plume,  its  latitude  and  longitude  of 
maximum  Al,  and  age  class.  Additional  back  trajectory  runs, 
using  starting  altitudes  of  0.5-  and  1-km  lower  than  the  maxi¬ 
mum  CALIPSO  plume  heights,  have  shown  that  the  age  class 
assigned  to  each  plume  is  robust. 
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Table  1.  The  dates  of  large-AI  biomass  burning  plumes  with  coin¬ 
cident  CALIPSO  data. 


Date 

Plume  height 
(km) 

Maximum 

Al 

Lat. 

Lon. 

Age 

class^ 

7/24/06 

4.5 

8.2 

62.5 

103.5 

Y 

7/25/06 

10.6 

9.2 

63.5 

100.5 

Y 

7/26/06 

12.0 

14.6 

61.5 

110.5 

Y 

7/27/06 

12.6 

11.1 

64.5 

114.5 

Y 

7/27/06 

7.5 

7.8 

72.5 

97.5 

Y 

9/05/06 

12.3 

11.4 

48.5 

-89.5 

Y 

9/05/06 

9.9 

8.6 

47.5 

-109.5 

Y 

12/19/06 

15.6 

9.3 

-37.5 

-78.5 

A 

12/19/06 

14.0 

7.6 

-39.5 

-69.5 

A 

5/26/07 

7.5 

8.1 

47.5 

-63.5 

Y 

6/10/07 

15.6 

15.4 

39.5 

122.5 

Y 

7/08/07 

12.0 

7.9 

33.5 

-104.5 

Y 

7/30/07 

8.7 

10.1 

49.5 

-103.5 

Y 

8/01/07 

6.0 

7.2 

47.5 

-94.5 

Y 

8/01/07 

4.5 

7.1 

-8.5 

11.5 

Y 

8/13/07 

4.8 

8.6 

47.5 

-109.5 

Y 

8/14/07 

7.5 

OO 

CO 

45.5 

-89.5 

Y 

8/15/07 

6.0 

7.0 

44.5 

-98.5 

Y 

8/17/07 

7.0 

7.2 

43.5 

-99.5 

Y 

7/02/08 

7.0 

7.3 

53.5 

147.5 

Y 

2/09/09 

13.8 

12.4 

-24.5 

-175.5 

Y 

2/11/09 

17.8 

7.0 

-32.5 

-169.5 

A 

2/12/09 

19.0 

7.1 

-30.5 

-175.5 

A 

“  Y  =  Young  (<2d  old)  and  A  =  Aged  (>2d  old) 


4  Results 

The  inclined  black  line  in  Fig.  3  represents  a  minimum  al¬ 
titude  of  the  top  of  a  plume  observed  with  a  given  AI  value 
(lower  bounding  line): 

Minimum  height  (km)  =  1 .2  AI  —  5.4  (1) 

This  bounding  line  suggests  that  for  all  plumes,  when  the  ob¬ 
served  AI  value  is  larger  than  9,  the  plume  altitude  is  higher 
than  5  km.  This  provides  an  AI  threshold  for  identifying 
high-altitude  injection  plumes  even  if  a  direct  height  mea¬ 
surement  is  not  available.  This  threshold  value  of  AI  is  sim¬ 
ilar  to  that  reported  by  Fromm  et  al.  (2008b),  who  reported 
that  unusually  large  AI  values  (10  or  larger)  were  associated 
with  smoke  in  the  UT/LS. 

While  a  threshold  AI  value  is  useful  for  categorizing  high 
altitude  plumes,  there  is  more  information  that  can  be  ex¬ 
tracted  from  Fig.  3.  The  red  line  is  the  best-fit  line  to  the 
young  plumes  with  a  corresponding  equation  of: 

Average  height  (km)  =  1 . 1  AI  —  1 .2  (2) 

There  is  a  clear  tendency  for  young  plumes  to  show  AI  val¬ 
ues  increasing  with  increasing  altitude.  The  correlation  co¬ 
efficient  between  AI  and  the  altitude  is  0.81.  This  suggests 
that  a  common  controlling  process  is  operative  in  all  young 
biomass  burning  plumes  identified  here.  Using  a  radiative 


transfer  model,  Wong  and  Li  (2002)  showed  that  the  top- 
of-atmosphere  reflectance  at  the  channel- 1  wavelength  (cen¬ 
tered  around  0.65  pm)  of  the  Advanced  Very  High  Resolution 
Radiometer  (AVHRR)  sensor  reached  a  constant  value  when 
aerosol  optical  depth  exceeded  a  threshold  value.  The  value 
of  aerosol  optical  thickness  at  which  the  asymptote  occurs 
seems  to  be  independent  of  the  value  of  SSA.  More  recent 
calculations  with  the  wavelengths  of  335.5  and  376.5  nm  (see 
Penning  de  Vries,  2010)  have  reached  a  very  similar  con¬ 
clusion.  We  hypothesize  that  most  of  the  analyzed  large-AI 
young  plumes  are  sufficiently  optically  thick  so  that  AI  is 
close  to  or  reaches  an  asymptote  with  respect  to  aerosol  op¬ 
tical  depth  at  ~360  nm,  similar  to  the  behavior  suggested  by 
Wong  and  Li’s  calculation  for  a  longer  wavelength  (2002). 
When  this  occurs,  and  because  the  SSA  is  not  likely  to  vary 
a  great  deal  (Jeong  and  Hsu,  2008),  the  plume  height  will  de¬ 
termine  the  value  of  AI.  That  is,  an  increasing  plume  height 
will  be  associated  with  a  monotonically  increasing  AI,  as 
shown  in  Fig.  3.  Figure  3  also  shows  the  variability  of  young- 
plume  heights  for  a  given  AI  value  (1 -sigma  R::  1.9  km).  This 
variability  may  reflect  physical  differences  between  plumes 
as  well  as  the  uncertainty  resulting  from  CALIPSO  measur¬ 
ing  only  a  small  portion  of  a  given  plume.  Uncertainty  in  the 
AI  measurements  themselves  to  variability  in  the  best-fit  re¬ 
lationship  is  problematic  to  assess  (De  Graaf  et  al.,  2005a,  b). 

The  four  aged  plumes  (four  or  five  days  old)  shown  in  blue 
in  Fig.  3  indicate  that  at  similar  altitudes,  aged  plumes  have 
lower  AI  values  than  young  plumes.  These  four  plumes  all 
originated  from  Australian  wildfires  (see  Table  1).  Mixing 
with  clean  air  as  the  plumes  aged,  was  likely  the  driving  fac¬ 
tor  in  reducing  the  optical  depth,  and  thus,  the  AI.  Previ¬ 
ous  studies  (Reid  et  al.,  1998;  Abel  et  al.,  2003)  also  reveal 
that  particle  aging  rapidly  increases  the  SSA  within  several 
hours  after  the  smoke  is  released  into  the  atmosphere.  Swap 
et  al.  (2003)  indicate  that  chemical  and  physical  evolution 
of  smoke  plumes  is  very  rapid  in  the  first  hour,  where  after 
properties  remain  relatively  stable  for  hours  to  days.  These 
studies  suggest  that  particle-aging  processes  may  have  a  mi¬ 
nor  role  in  reducing  the  AI  value  for  plumes  older  than  sev¬ 
eral  hours.  Two  of  the  aged  plumes  (with  heights  of  14  and 
15.6  km)  were  from  a  fire  on  14  December  2006  (Dirksen  et 
al.,  2009),  while  the  other  two  plumes  (17.8  and  19  km)  were 
from  a  fire  on  7  February  2009.  These  aged  plumes  were 
visible  in  the  AI  data  for  several  days.  Therefore  we  were 
able  to  record  their  AI  values  on  their  first  observation  day 
and  estimate  the  corresponding  heights  using  the  trajectory 
model.  The  model  was  run  forward  starting  from  15  De¬ 
cember  2006  and  8  February  2009  at  day-1  plume  locations. 
The  trajectories  starting  at  13.5  and  16.5  km  pass  through  the 
CALIPSO-observed  aged-plume  locations  and  top  heights, 
indicating  the  plume  tops  are  most  likely  located  at  13.5  km 
for  the  2006  case  and  16.5  km  for  the  2009  case  when  they 
were  young  (1  day  old).  The  inferred  young-plume  data  (red 
triangles)  were  added  into  the  Al-plume-height  plot  as  shown 
in  Fig.  4.  The  four  aged  plumes  (blue)  were  inside  or  near 
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Fig.  4.  Scatter-plot  of  OMI  AI  and  CALIPSO  maximum  plume 
height.  Red  and  blue  symbols  represent  “young”  and  “aged” 
plumes,  respectively,  as  described  in  Fig.  3.  Solid  shaded  area 
shows  ±2  standard  deviation  (2-sigma  3.8  km)  across  the  best  fit 
of  the  young  plumes  indicated  by  red  crosses.  The  young-plume 
data  in  red  triangles  were  inferred  from  forward  trajectory  model 
calculations  and  OMI  observations  3  or  4  days  before  the  obser¬ 
vations  (blue  crosses).  The  four  aged  plumes  (4-5  days  old)  were 
inside  or  near  the  dash  shaded  area  (2-sigma)  when  they  were  young 
(1  day  old). 

the  2-sigma  shaded  area  and  had  a  lower  altitude  when  they 
were  young  (1  day  old).  This  indirectly  implies  that  dilution 
processes  resulting  in  much  smaller  optical  depths  are  one  of 
the  causes  for  aged  plumes  deviating  from  the  Al-height  re¬ 
lationship  of  young  plumes.  This  also  indirectly  confirms  the 
findings  of  Real  et  al.  (2007)  who  showed  the  important  role 
of  dilution/mixing  processes  on  reducing  simulated  pollutant 
levels  (or  concentrations). 

5  Finding  the  geographical  distrihution  and  historical 
number  of  high-altitude  plumes 

The  knowledge  of  the  geographical  distribution  and  histori¬ 
cal  number  of  high-altitude  biomass  burning  plumes  is  criti¬ 
cal  for  further  understanding  of  their  associated  radiative  and 
climate  effects.  In  this  section,  we  use  the  threshold  value  of 
AI>9  and  lower  bound  Eq.  (1)  to  identify  plumes  above  5  km 
and  find  the  geographic  distribution  and  historical  number  of 
high-altitude  plumes. 

The  global  distribution  of  all  AI>9  plumes  during  the  en¬ 
tire  data  record  is  displayed  in  Fig.  5.  Given  the  lower  bound¬ 
ing  line  shown  in  Fig.  3,  these  plumes  are  expected  to  be  at 
altitudes  of  5  km  or  higher  and  represent  the  total  number 
of  observations.  That  is,  these  plumes  have  been  counted 
each  day  they  were  observed.  This  total  number  is  an  essen¬ 
tial  parameter  when  estimating  the  influence  of  pollutants  on 


the  UT/FS.  Sources  of  error  include  the  omission  of  high- 
altitude  plumes  due  to  optical  thinness  and  the  limited  tem¬ 
poral  resolution  of  polar-orbiting  satellites.  A  careful  dust¬ 
screening  process  has  helped  reduce  the  possibility  of  con¬ 
tamination  by  dust  plumes. 

North  America,  Russia/Northeast  Asia,  and  Australia  are 
three  major  biomass  burning  source  regions  yielding  large 
AI  values.  The  number  of  large  AI  plumes  from  each  of 
these  three  regions  during  the  entire  AI  record  is  shown  in 
Fig.  6.  The  number  of  the  dust  plumes  screened  out  for  Rus¬ 
sia/Northeast  Asia  is  also  displayed  in  the  figure.  A  few  re¬ 
moved  cases  which  possibly  included  smoke  mingled  with 
dust  do  not  affect  the  following  conclusion.  Typically,  large 
AI  plumes  occur  more  frequently  over  North  America  (75) 
than  over  Russia/Northeas  Asia  (60)  or  Australia  (28).  Pre¬ 
vious  studies  have  suggested  that  Russian  fires  are  largely 
lower  intensity  fires,  burning  mainly  surface  fuels  (Furyaev, 
1996;  Conard  and  Ivanova,  1997;  Wooster  and  Zhang,  2004), 
whereas  in  North  America  much  higher  intensity  “crown” 
fires  predominate  (Conny  and  Slater,  2002;  Wooster  and 
Zhang,  2004).  Crown  fires  have  a  larger  potential  to  inject 
pollutants  to  higher  altitudes,  which  is  consistent  with  our 
current  analysis.  Although  originating  from  a  much  smaller 
forest  area,  the  number  of  Australian  high-altitude  plumes 
is  also  noteworthy.  These  occurrences  mainly  resulted  from 
the  bushfires  in  Victoria,  Australia  during  December  2006 
and  around  7  February  2009.  An  increase  in  fire  number,  du¬ 
ration,  and  intensity  in  Australia  has  been  predicted  by  pre¬ 
vious  researchers  (Beer  and  Williams,  1995;  Williams  et  al., 
2001;  Cary,  2002).  It  will  be  interesting  to  see  in  the  coming 
years  if  this  trend  persists. 

Although  Africa  is  the  single  largest  continental  source  of 
biomass  burning  emissions  (Robert  et  al.,  2009),  large  AI 
plumes  are  absent  in  the  entire  climatology.  Fires  in  Africa 
burn  biomass  such  as  savanna  and  grassland,  which  produce 
broader  and  lower  injection  layers  due  to  lower  fire  inten¬ 
sity  (Freitas  et  al.,  2007).  These  lower  plume  layers  result  in 
prevalent  small-AI  plumes  over  Africa. 

6  Conclusions 

This  study  shows  the  AI  data  set  can  be  used  to  identify  high- 
altitude  fire  plumes.  The  derived  best-fit  relationship  (Eq.  2) 
between  AI  and  maximum  plume  height  for  young  plumes 
(within  2  days  of  the  fire)  can  help  validate  the  vertical  place¬ 
ment  of  smoke  plumes  in  chemical  transport  models  and  can 
also  provide  a  method  for  the  near-real-time  monitoring  of 
high-altitude  events. 

For  fire  plumes  measured  with  coincident  OMI  and 
CAFIPSO  data  from  June  2006-February  2009,  we  found 
the  plume  height  is  >5  km  when  the  AI  value  is  >9. 
During  the  period  1978-2009,  large  AI  (>9),  and  there¬ 
fore  likely  high-altitude,  plumes  occurred  more  frequently 
over  North  America  (~75)  than  over  Australia  (~28)  or 


Atmos.  Chem.  Phys.,  10,  6461-6469,  2010 


WWW.  atmos-chem-phy  s  .net/ 1 0/646 1/2010/ 


H.  Guan  et  al.:  A  multi-decadal  history  of  biomass  burning  plume  heights 


6467 


Fig.  5.  Global  distribution  of  high-altitude  plumes  with  an  AI  value  >9  (i.e.  height  >5  km)  during  the  entire  AI  data  record  (31  October  1978 
to  13  February  2009).  181  plumes  were  identified,  and  163  of  those  are  located  in  these  three  main  regions. 


120 

100 


Russia  and  Northeast  Asia  North  America  Australia 

Fig.  6.  Number  of  large-AI  (>9)  observations  over  the  three  regions 
(shaded  hars)  shown  in  Fig.  5  during  the  entire  AI  data  record.  The 
white  har  represents  the  number  of  dust  plumes  screened  out  over 
Russia/Northeast  Asia. 


Russia/Northeast  Asia  (~60).  According  to  the  derived  re¬ 
lationship,  during  this  time  interval,  181  plumes,  in  vari¬ 
ous  stages  of  their  evolution,  reached  altitudes  above  8  km. 
Considering  their  large  mass  loading  and  long  lifetime,  these 
high-altitude  plumes  may  have  the  most  important  impact  on 
the  radiative  budget,  climate,  and  chemistry.  Further  tracing 
their  dilution  and  aging  using  satellite  data  may  give  addi¬ 
tional  information  regarding  their  transport  and  fate. 

This  large  AI  screening  approach  provides  a  useful  and 
fast  tool  to  estimate  the  approximate  top  height  of  plumes. 
However,  this  height  determination  approach  is  biased  in  fa¬ 
vor  of  dense  high-altitude  smoke  because  of  the  reduced  sen¬ 
sitivity  of  AI  to  low-altitude  aerosols.  Methods  for  detecting 
lower  plumes  and  optically-thin  high-altitude  plumes  are  be¬ 


ing  explored.  The  effectiveness  of  the  linear  relationship  be¬ 
tween  AI  and  maximum  plume  height  for  multi-layer  plumes 
or  for  thick  (a  few  km)  plumes  with  a  large  difference  be¬ 
tween  the  maximum  and  effective  heights  also  needs  to  be 
clarified  in  the  future. 
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